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ABSTRACT 

Context. Infrared-faint radio sources (IFRS) form a new class of galaxies characterised by radio flux densities between tenths and tens 
of mjy and faint or absent infrared counterparts. It has been suggested that these objects are radio-loud active galactic nuclei (AGNs) 
at significant redshifts (z > 2). 

Aims. Whereas the high redshifts of IFRS have been recently confirmed based on spectroscopic data, the evidence for the presence of 
AGNs in IFRS is mainly indirect. So far, only two AGNs have been unquestionably confirmed in IFRS based on very long baseline 
interferometry (VLBl) observations. In this work, we test the hypothesis that IFRS contain AGNs in a large sample of sources using 
VLBl. 

Methods. We observed 57 IFRS with the Very Long Baseline Array (VLBA) down to a detection sensitivity in the sub-mJy regime 
and detected compact cores in 35 sources. 

Results. Our VLBA detections increase the number of VLBl-detected IFRS from 2 to 37 and provide strong evidence that most—if 
not all—IFRS contain AGNs. We find that IFRS have a marginally higher VLBl detection fraction than randomly selected sources 
with mJy flux densities at arcsec-scales. Moreover, our data provide a positive correlation between compactness—defined as the 
ratio of milliarcsec- to arcsec-scale flux density—and redshift for IFRS, but suggest a decreasing mean compactness with increasing 
arcsec-scale radio flux density. Based on these findings, we suggest that IFRS tend to contain young AGNs whose jets have not formed 
yet or have not expanded, equivalent to very compact objects. We found two IFRS that are resolved into two components. The two 
components are spatially separated by a few hundred milliarcseconds in both cases. They might be components of one AGN, a binary 
black hole, or the result of gravitational lensing. 

Key words. Techniques: interferometric - Galaxies: active - Galaxies: high-redshift - Galaxies: nuclei - Radio continuum: galaxies 


1. Introduction 

One of the most controversially discussed topics in current as¬ 
trophysics is the question of how active galactic nuclei (AGNs) 
influence star forming activity in their host galaxies. This in¬ 
teraction is known as feedback, and both negative and positive 
feedback of the AGN have been proposed. 

Negative AGN feedback is thought to be caused by the AGN 
heating or disrupting the surrounding gas and, by this, quenching 
star formation. This scenario consists of two consecutive phases 
and may be caused by a merger of two galaxies. In the first phase, 
which is called cold mode (or quasar mode; e.g. Kauffmann et al. 
|2003[[&oton et al.|2006[[Hardcastle et al.|20()7| l, cold gas from 

the galaxies involved in the merger fuels the accretion accom¬ 
panied by high star forming activity. Finally, the AGN heats or 
blows away the remaining gas and star formation ends. This rep¬ 
resents the transition to the second phase, the so-called hot mode 
accretion (also known as radio mode; e.g. |Narayan & Yi|1995| l. 
In this phase, the supermassive black hole (SMBH) can only ac¬ 


crete hot gas from the halo in so-called advection dominated ac¬ 
cretion flows. 

The accretion in the cold mode phase is radiatively efficient 
and observational evidence is difficult because of the highly ob¬ 
scured nucleus. In contrast, accretion in the hot mode phase is 
typically well below the Eddington limit and the accretion flow 
radiatively inefficient. In hot mode, AGN feedback is kinetic 
and the total energy output usually dominated by the mechan¬ 
ical power of the radio jets. Radio jets are associated with both 
accretion modes and the radio emission can be used to trace the 
mechanical jet power (e.g. [Alexander & Hickox|2012| and refer¬ 
ences therein). 

This scenario of negative AGN feedback is an important in¬ 
gredient for the current preferred cosmological ACDM model. 
In this model, negative AGN feedback is needed to make the 
number of massive and luminous galaxies in simulations consis¬ 
tent with the observed number (|Springel et al.|2005[ [Croton et al.| 
[2006l l. 

Contrary to this widely accepted negative AGN feedback, 
examples have been found where AGN activity enhances star 
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formation, i.e. showing a positive feedback (e.g. Klamer et al.| 
|2004|[Gaibler et al.|2012[|Zinn et al.|2013||Karouzos et al.|2014| ^ 

These observations are explained by the AGN jets propagating 
through the gas of the host galaxy, generating shocks that trig¬ 
ger the gravitational collaps of the gas and, by this, boost star 
formation. 

One important class of object in studying AGN feedback 
processes is the class of high-redshift radio galaxies (HzRGs) 
since these galaxies can be observed out to high redshifts. 
HzRGs are very powerful radio-galaxies that contain signihcant 
emission both from a starburst activity and AGN activity. They 
are expected to be the progenitors of the most massive galaxies 
in the local universe (e.g. [Seymour et al.|2007[|De Breuck et al.| 
|2010| l. Only around 200 of these objects—which are defined by 
z > 1 and L 3 GHZ > lO^^WHz '—are known in the entire sky. 
However, a new class of object has recently been found that sug¬ 
gests a link to HzRGs: the class of infrared-faint radio sources. 

1.1. Discovery and definition of infrared-faint radio sources 

Infrared-faint radio sources (IFRS) are characterised by radio 
emission of the order of tenths to tens of mJy and associated deep 
near-infrared faintness. [Norris et al.| ( |2006| ) and [Middelberg et ah] 
( 2008a[ l discovered these objects in the deep radio maps of the 
Australia Telescope Large Area Survey (ATLAS) in the Chandra 
Deep Field South (CDFS) and the European Large Area IR space 
observatory Survey South 1 (ELAIS-Sl) as lacking infrared (IR) 
counterparts in the co-located Spitzer Wide-area Infrared Ex- 
tragalactic Survey (SWIRE; [Lonsdale et al.[ 2003 [ l. [Zinn et 3T] 
([201 l|l dehned two survey-independent criteria for the selection 
oflERS: 

(i) s iAGHzlS3.6pm > 500 and 

(ii) < 30//Jy. 

The high radio-to-IR flux density ratios, ensured by the hrst cri¬ 
terion, show that lERS are clear outliers from the radio-to-IR 
correlation. The second criterion is equivalent to a distance se¬ 
lection and prevents ordinary objects at z < 1.4 from being in¬ 
cluded in this class of object. 

Based on the two selection criteria, [Zinn et al^pOl l[ l com¬ 
piled a catalogue of 55 lERS in the deep helds of the CDES, 
ELAIS-Sl, Spitzer extragalactic Eirst Look Survey (xELS), and 
the Cosmological Evolution Survey (COSMOS), based on the 
work by Norris et ari ( |2006| l, [Middelberg et al. ( [2008a| l, and [Garn[ 
[& Alexander[ ( [20()8| . Later, lERS were also found in the Euro¬ 
pean Large Area IR space observatory Survey North 1 (ELAIS- 
N1) held by [Banheld et al.[ ( |201 l[ l and in the Lockman Hole held 
by Maini et al. (submitted). Around 100 lERS have been found 
in these deep helds, covering a total area of around 35 deg^. 

Recently, [Collier et al. ( [2014 1 , for the hrst time, used a differ¬ 
ent approach and looked for lERS in much shallower radio and 
IR data, which covered a much larger area compared to the deep 
helds mentioned above. ICollier et al.l used data from the Uni- 
hed Radio Catalogue (URC; [Kimball & Ivezic][2008[ [Kimball[ 
[& Ivezic|2014 1 and from the all-sky Wide-Eield IR Survey Ex¬ 
plorer (WISE; Wright et al.|2010 i. Based on these data, [Collier[ 
[et al.[ compiled a catalogue of 1317 lERS, all of them fulhlling 
both selection criteria from|Zinn et al.[([201 l|l. 


1.2. The properties of IFRS 


shifts (z > 2). Different studies found evidence for this sugges¬ 
tion. Norris et akl l 2007[ l and Middelberg et al. (2008b I presented 
very long baseline interferometry (VLBI) observations of lERS 
and detected compact cores in two lERS. [Norris et alij used the 
Australian Long Baseline Array (LBA) and targeted two lERS 
from the CDES, out of which one was detected, showing a total 
VLBI Hux density of 5.0 mJy at 1.6 GHz. The source was un¬ 
resolved on VLBI scales and its size measured to be less than 
0.03", corresponding to less than 260 pc at any redshift. Mid 


[delberg et al.[ also used the LBA and observed four lERS from 
the ELAIS-Sl held at 1.6GHz or 1.4GHz. One out of these 
four lERS was detected with a VLBI Hux density of 12.5 mJy at 
1.6 GHz. Based on a Hux density of 7 mJy on the longest base¬ 
lines, Middelberg et al. [ measured a lower limit on the brightness 
temperature of 3.6 x 10® K. Three lERS were included in the 
sample observed with the Very Long Baseline Array (VLBA) in 
the CDES by Middelberg et al.| ( |201 la[ l; however, they remained 
undetected. The two VLBI detections of lERS showed that at 
least a fraction of lERS contain AGNs. 


Garn & Alexander (2008 1 and Huynh et al. (2010 1 found that 


obscured star forming galaxies cannot reproduce the character¬ 
istics of lERS because lERS clearly deviate from the radio-IR 
correlation, providing another hint for the AGN content in lERS. 
Eurthermore, [Garn & Alexander[ and [Huynh et al. [ presented the 
hrst spectral energy distribution (SED) modelling of lERS and 
showed that 3C sources like 3C 273 are in agreement with the 
characteristics of lERS if these sources are at redshifts z > 2. 
[Herzog et al.[ ( [2014[ ) showed that the SED of their sample of 
lERS can only be explained by radio-loud AGN templates. 

A potential link between lERS and HzRGs was hrst sug¬ 
gested by Huynh et al. ( 2010|l based on the sim ilarly high radio- 
to-IR hux density ratios. [Middelberg_ et al.[ ( [201 Ib[ ) showed 
that lERS have steeper radio spectr^ (median spectral index 
a = -1.4) than the general radio source population (a - -0.86) 
and the AGN source population (a - -0.82), using data in the 
ELAIS-Sl held between 2.3 GHz and 8.4GHz. Moreover, [Mid-[ 
delberg et al.[ found that the radio spectrum of lERS is even 
steeper than that of HzRGs (a - -1.02). [Norris et aL] ( [201l| l 
pointed out that HzRGs are the only objects known at high red¬ 
shifts that share the extreme radio-to-IR Hux density ratios with 
lERS. Based on this similarity and the deep IR faintness of lERS, 
[Norris et al.[ suggested that lERS might follow the correlation 
between 3.6 pm hux density and redshift found for HzRGs ([Sey¬ 
mour et ar][2007|l, similar to the K - z relation ([Willott et al. 


2003 1 . Herzog at al. (submitted) showed that the non-detection 


of lERS in deep far-IR Herschel observations can only be ex¬ 
plained by SED templates of HzRGs. 

Collier et al. (20141 and Herzog et al. (2014|l presented the 


Since the hrst detection of lERS by Norris et al. ( [2006 j l it has 
been suggested that lERS are radio-loud AGNs at signihcant red- 


hrst spectroscopic redshifts of lERS and found all of 22 but 
one—which is most likely a misidentihcation or a star forming 
galaxy with an AGN—redshifts in the range 1.7 < z < 3.0, con- 
hrming the suggested high-redshift nature of lERS. Eurthermore, 
both studies found their data in agreement with the suggested 
correlation between near-IR hux density and redshift, indicat¬ 
ing that most lERS in deep helds—which were summarised by 
[Zinn et aL| ( [201 ip and Maini et al. (submitted)—might be at even 
higher redshifts. 

1.3. Populations of IFRS found in deep and shallow surveys 

Here, we consider the relationship between lERS found in deep 
surveys of small area and those found in shallow all-sky surveys. 


The spectral index is dehned as 5 oc v“ 
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First, all IFRS found in the various works mentioned above fulfil 


the—in some cases slightly changed—selection criteria by Zinn 
et al.|(20TT]) and therefore qualify as IFRS. Maini et al. (submit 


ted) lowered the radio-to-IR flux density criterion and replaced 
the IR flux density criterion by an extension criterion. Others— 
like Norris et al. ( |2007 1 and Middelberg et al. (2008a i—required 
IFRS to be undetected at 3.6/tm without applying a radio-to-IR 
flux density criterion. Nevertheless, these IFRS, which were se¬ 
lected based on slightly different criteria, are close to fulfil the 
criteria by |Zinn et al.| 

While the IFRS in deep fields usually have 1.4 GHz flux den¬ 
sities of tenths of mJy to a few mJy, the IFRS in the all-sky sur¬ 
vey have median radio flux densities of several tens of mJy, some 
reaching even several hundred mJy. Similarly, each IFRS in the 
sample from |Collier et al.| ( |2014) l provides an IR counterpart at 
3.4yum with a mean flux density of around 25 pJy, while a signif¬ 
icant fraction of IFRS in deep fields has no IR counterpart in the 
even deeper 3.6/rm data. Thus, the median IFRS in deep fields is 
both radio and IR-fainter than the median IFRS in the shallower 
all-sky sample from |Collier et al.| 

All spectroscopic redshifts of IFRS from both the deep fields 
and the shallow survey were found to be in the same redshift 
range 1.7 < z < 3.0. However, a selection effect putatively biases 
the observed spectroscopic redshift distribution of IFRS in the 
deep fields since the IFRS with known spectroscopic redshifts 
in these fields are the optically and IR brightest IFRS in that 
sample. Since their optical, IR, and radio properties are similar to 
those of the all-sky IFRS, it is expected that these IFRS represent 
the overlap between the fainter IFRS population in deep fields 
and the brighter all-sky IFRS population (Herzog et al.|2014)l. 


It has been suggested by Collier et al. (|2014|l that their 


IFRS sample consists of the lowest-redshift IFRS while the IFRS 
found in deep fields are on average at higher redshifts. This sug¬ 
gestion is in agreement with the correlation between 3.6yum flux 
density and redshift discussed by |Norris et al.| (|201 1||, |Collier| 
|et al.| and [Herzog et al.| ( |2014| l. The overlapping spectroscopic 
redshifts found for both subsets are in agreement with this sug¬ 
gested unification, too. 

Collier et al. (20141 could only set a lower limit of ~ 
0.1 deg“ for the sky density of IFRS with S i. 4 ghz ^ 7.5 mJy 
because of the non-uniform sensitivity of the WISE survey. In 
contrast, the sky density of IFRS in deep fields is o f the order 
of a few per square degree and might reach 30 deg“^ (Zinn et al. 

|20TT) . 

In this paper, we test the AGN content in IFRS based on 
VLBI observations with the VLBA of a large number of sources 
taken from the all-sky catalogue of IFRS ( [Collier et al.|2014 1. In 
Sect. 1^ we describe our sample and the observing strategy. We 
discuss data calibration, imaging, flux measuring, and redshifts 
in Sect.[^ We analyse our data in Sect. with respect to detec¬ 
tion fraction (Sect. pTTj ), compactness (Sect. [4~2] l, and individual 
sources (Sect. [4.3] l. We discuss the implications of our analysis in 
Sect. and present our conclusions in Sect. Throughout this 
paper, we use flat ACDM cosmological parameters Qa = 0.7, 
Qm = 0.3, Ho — 70 km s * Mpc ' and the calculator by Wright 
(2006[l. We quote 1 cr confidence intervals of binomial population 


proportions based on the bayesian approach, following|Cameron| 

( [ 20 ^ . 


2. Sample and observations 


We selected all IFRS from the catalogue from Collier et al. 


([2014[l which were located within 1 deg of a VLBA calibrator. 


This low angular separation between source and calibrator en¬ 
sured the phase coherence required for VLBI observations. Since 
1.4 GHz VLBI observations of the calibrators were not available, 
we required the calibrators to have a 2.3 GHz flux density of at 
least 0.2 Jy on a baseline of 5000 km. Out of the 1317 ILRS pre¬ 
sented by Collier et al. 110 were found to provide a calibrator 


which fulfills the given conditions. 


A VLBI detection provides unambiguous evidence for an 
AGN because compact radio emission in AGNs is a non-thermal 
process and results in brightness temperatures of more than 
10® K to which our observations are sensitive. In contrast, com¬ 
pact radio emission in starburst galaxies, which is usually dom¬ 
inated by thermal free-free emission, i s represented by br ight- 
ness temperatures of around 10"^ K (e.g. 


Condon et al. 


19911. Al¬ 


though brightness temperatures of 10® K can also be produced 


by very luminous radio supernovae (SNe; [Huang et al. 1994[ 


[Smith et al.|1998[ ), [Kewley et al.[ ( [2000| l showed that the proba 
bility for a VLBI-detected radio SN in a galaxy sample is very 
low. Therefore, a VLBI detection provides strong evidence for 
an AGN. However, it should be noted that the reverse is not true, 
i.e. the non-detection of a source in a VLBI observation does 
not imply the non-existence of an AGN. Instead, a VLBI non¬ 
detection implies significant extended emission compared to the 
compact core. The ratio of extended emission to core emission 
depends on beaming which can boost or suppress the compact 
core emission, AGN age, and the surrounding medium, affecting 
the brightness and extent of the diffuse radio lobes. 


Since this is a detection experiment, the wv-coverage is not 
critically important and a minimum number of six out of ten 
VLBA antennas was requested. Since the individual observa¬ 
tions were short and independent of the weather conditions be¬ 
cause of the observing frequency of 1.4 GHz, the observations 
were scheduled in filler time. Although the maps resulting from 
the data will be of rather poor quality, they will unambiguously 
resolve potential compact components. 


Out of 110 proposed objects, 57 ILRS were observed in 
semester 14A in project BH197. The 57 observed ILRS were 
randomly selected based on the ILRS positions and available 
filler time at the VLBA. These observed sources have 1.4 GHz 
integrated flux densities between llmJy and 183 mJy in the 
NRAO VLA Sky Survey (NVSS; [Condon et ah][19^ . The 
VLBA observations were set to a bandwidth of 32 MHz in each 
of the eight basebands which were observed in dual polarisation 
at 1.4 GHz, resulting in a total data rate of 2048 Mbps. 


Each of the 57 epochs had a total observing time of one hour. 
We decided to use two different approaches for the scan settings, 
depending on the distance between source and calibrator. If the 
separation between source and calibrator was less than 25', we 
continuously pointed at the position in between ILRS and cali¬ 
brator to prevent unnecessary nodding between the two sources. 
If the separation was more than 25', we alternately observed the 
calibrator for 60 s and the source for 225 s, starting and ending 
with a scan on the calibrator. The resulting observing time on 
the ILRS was around 45 min. The data were correlated using the 
VLBA Distributed LX (DiLX) software correlator ([Deller et al.| 
[20071 [20 TT] i. 
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3. VLBA data calibration, optical properties, and 
redshifts 

3.1. Data calibration, imaging, and flux measurement of the 
VLBA data 

Calibration and imaging of the individual epochs was carried out 


based on a ParselTongue script. ParselTongue (Kettenis et al. 
|2006| l is a Python-based interface to the Astronomical Image 
and Processing Systenj^ (AIPS). The calibration and imaging 
strategy used here is very similar to the procedure described by 
|Deller & Middelberg| ( |2014] l. However, we calibrated the ampli¬ 
tudes using the technique suggested by VLBA Scientific Memo 
#37 (R. Craig Walker; Dec 15, 2014). 

In the pipeline, we loaded the data into AIPS and applied a 
priori flags and manual flags. We then corrected for ionospheric 
effects using the task TECOR, applied the latest earth orienta¬ 
tion parameters, and corrected for parallactic angles, in the latter 
cases using the task CLCOR. Amplitudes were calibrated based 
on autocorrelation data using the task ACCOR. We corrected for 
primary beam effects using the task CLVLB, following the pro¬ 


cedure outlined by Middelberg et al. (20131. Delay correction 


was carried out based on the calibrator, using the task FRING 
and a solution interval of 2 min. We applied a bandpass correc¬ 
tion using the task BPASS. The new task ACSCL, implemented in 
AIPS in consequence of the flux density calibration errors de¬ 
scribed in the VLBA Scientific Memo #37, kept the calibrated 
autocorrelation values at unity. The flux density calibration was 
completed by using APCAL, calibrating the amplitudes based on 
system temperatures and gains. We used the AIPS/ParselTongue 
implementation of Pieflag ( |Middelberg|2006 1 , dynspec-flagger, 
to automatically flag data affected by interference. Using the 
task CALIB, we performed one iteration of phase and ampli¬ 
tude self-calibration on the calibrator and applied the solution 
on the target. Finally, we imaged the UV data based on the 
task IMAGR into a map of 2048x2048 pixels with a pixel size 
of 1 milliarcsecond (mas)—matching the angular resolution of 
~ 5 mas—using uniform weighting. In this step, we cleaned the 
dirty image in a given box, down to a flux of three times the rms 
or up to 500 iterations. 

The mean synthesised beam size in our observations is 14.3 x 
4.7 mas^. Because the linear scale is limited in a ACDM cosmol¬ 
ogy at redshifts 0.5 < z < 12 between 4kpc/"and 8.5kpc/", we 
were able to convert this beam size over this redshift range to an 
area of between 57 pc X 19pc and 122pcx40pc. Therms in the 
final maps is around 60//Jybeam ', depending on the number of 
antennas. Of the 57 observations, 31, 21, 4, and 1 were carried 
out with 10, 9, 8, and 7 antennas, respectively. Since the longest 
VLBA baseline of 8611 km—^between the stations Mauna Kea 
and St. Croix—was available in all observations, the east-west 
angular resolution of around 5 mas is similar for all 57 observa¬ 
tions, whereas the north-south angular resolution varies slightly 
because of changing antenna availability. As an example. Fig. 
shows the flnal map of IFRS F0398. 

We measured the flux density of each component using the 
AIPS task 3MFIT, setting a box of size of 40 pixels around the 
components. For a few sources, image artefacts resulted in un¬ 
reasonable fluxes. In these cases, we manually measured flux 
densities using the task TVSTAT. Checks on control sources re¬ 
sulted in consistent flux measurements based on TVSTAT and 
JMFIT. For the brighter component in IFRS F0030, we measured 


the flux density using the blobcat package ( Hales et al.|2012|l 
because of the complex structure of this source. Following|Deller| 


l^http : //WWW. aips. nrao. edu/ 
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Fig. 1. VLBA map of IFRS F0398. The source is slightly resolved with 
a peak flux density of 13.7 mjy bearn^* and an integrated flux density of 
17.2 mJy. 


& Middelberg (2014|l, we set a 6.15cr detection limit for all 


sources, corresponding to a mean detection sensitivity of around 
450/tJy beam^'. Most of our sources are slightly resolved, but 
insufficiently resolved to determine the morphology. The result¬ 
ing flux densities und flux density upper limits are summarised 
in Table Since observations and data calibration are very sim¬ 
ilar to the approach taken by Deller & Middelberg who targeted 
more than 20 000 sources in their mJy Imaging VLBA Explo¬ 
ration (mJIVE) survey, we expect our flux densities to be of sim¬ 
ilar accuracy of 20%. 


3.2. Optical properties and redshifts 


Our sample was taken from the all-sky lERS catalogue from 
ICollier et al.|P014l l who cross-matched their sources with the 
Sloan Digital Sky Survey (SDSS) DR9 ( |Ahn et ^|2012| l. Eor 
the purpose of this work, we cross-matched our sample of 
VLBA-observed lERS to the recent SDSS DRIO ( |Ahn et ah] 


201411. 53 (93+;%) out of our lERS are covered by SDSS DRIO. 


Eleven (2U;%) out of these 53 lERS provide photometric coun¬ 
terparts which are all close to the sensitivity limit of SDSS. 


We used the software EAZY ( [Brammer et al.j [2008) 1 with 
the standard template set to derive photometric redshifts for our 
sources. Since most lERS in the catalogue from 
( |2014) l are only detected in the WISE bands W1 
W2 (4.6 pm )—apart from the radio detections which are not used 
in the redshift fitting—measuring photometric redshifts is im¬ 
possible for most of the lERS in our sample. However, mea¬ 
suring photometric redshifts is possible for those lERS with 
SDSS counterparts. Eor these eleven lERS in our sample with 
SDSS counterparts, we obtained ten photometric redshifts using 
EAZY. As examples, we show the resulting fits for lERS E0197 
and E0273 in Eig.j^ The fitting of lERS E0277 failed. Table 
summarises the photometric redshifts and—where applicable— 
spectroscopic redshifts for our subsample of lERS with SDSS 
DRIO counterparts. Eig.[^ shows the photometric redshifts as a 
function of the spectroscopic redshifts for the three lERS in our 
sample for which SDSS DRIO provides spectroscopic redshifts. 


Collier et al. 
(3.4 um) and 


Based on the SEDs resulting from the fitting of photomet¬ 
ric redshifts, we And that most of our fitted lERS are very blue. 
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Fig. 2. Resulting photometric redshift fit from EAZY for IFRS F0197 (left) and F0273 (right). The SED (orange line) shows the best fit template 
based on the peak redshift in the observer’s frame. Black dots with error bars represent photometric data points of the IFRS, whereas black 
triangles show the flux density of the fitted SED at the same wavelengths. The subplots in the upper right of both plots show the redshift-probability 
distribution. The redshift of the peak in the probability distribution is quoted below each subplot. 


Table 2. Redshift information for those lERS with SDSS DRIO detec¬ 
tions. Listed is the lERS ID, the u band model magnitude from SDSS 
DRIO, the photometric redshift measured in this work using EAZY, and 
the spectroscopic redshift from SDSS DRIO. The SED fitting for F0277 
failed. 


IFRS 

u 

■^phot 

■^spec 

ID 

[mag] 



F0146 

20.21 + 0.05 


- 

F0194 

23.82 + 0.95 

0.78!°;« 

- 

F0197 

22.83 + 0.32 

o C4+0.14 

2.1150 + 0.0014 

F0273 

20.84 + 0.08 

1 65 

- 

F0277 

24.16 + 0.86 

- 

- 

F0293 

22.94 + 0.33 

o (AO+0.09 
■^■^^-0.08 

- 

F0398 

20.49 + 0.06 

o ^4+0.07 
^■^^-0.08 

2.55265 + 0.00021 

F0726 

23.21+0.52 

3 Q7+0.26 
-0.24 

- 

F0732 

22.88 + 0.27 

2 26 

^■^^-0.28 

- 

F0912 

21.55 + 0.10 

o 47+0.06 
-0.06 

2.61873 + 0.00023 

F1037 

25.14 + 0.79 

^■*^^-0.08 

- 


These characteristics can be explained by a type IAGN in the op¬ 
tical. However, we stress that we might be significantly affected 
by selection bias since these IFRS are very close to the detection 
sensitivity of SDSS and SDSS is more sensitive to blue objects. 
Nevertheless, this study shows that at least some IFRS are very 
blue, non-dusty galaxies. 


4. Analysis 

4. 1. VLB! detection fraction 

In our VLBA observations, we detected 35 (61^®%) out 
of 57 observed IFRS, showing peak flux densities between 
0.4mJybeam“' and 26.4mJybeam“' and integrated flux densi¬ 
ties between 0.6 mJy and 40.4 mJy as listed in Table 

The detection of a source in VLBI observations with bright¬ 
ness temperatures above 10® K—which is reached for our VLBA 
observations—is an unambiguous sign for an AGN as discussed 



Fig. 3. Comparison of the photometric redshifts obtained from EAZY 
and the spectroscopic redshifts from SDSS DRIO. The error bars show 
the Itr uncertainties in the photometric redshifts as determined by 
EAZY. The solid line represents the positions of exact agreement be¬ 
tween photometric and spectroscopic redshifts. 


in Sect.|^ However, we note that the reverse is not true, i.e. the 
non-detection of a source in VLBI observations does not exclude 
the existence of an AGN. The detection of 35 out of 57 IFRS in 
our VLBA observations provides strong evidence that most—if 
not all—IFRS contain AGNs. 


4.1.1. VLBI detection fraction compared to other samples 

The detection fraction in our VLBA observations of 61;^®% down 
to a 6.75cr detection limit of ~ 0.45 mJy beam^' is significantly 
higher than the detection fractions found by [Garrington et al.| 
( | 1999) 1 or |Deller & Middelberg ( 2014| l who targeted large sam¬ 
ples of sources from the Faint Images of the Radio Sky at Twenty 
centimetres (FIRST; |Becker et al.|1995| l survey in VLBI obser- 
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vations and detected 35% and 20%, respectively. However, the 
sensitivity of the respective VLBI observations and the sample 
selection criteria were different. [Garrington et al. had a detec¬ 
tion sensitivity between 1 mJy and 2 mJy and targeted FIRST 
sources with 1.4 GHz peak flux densities above lOmJy, whereas 


Deller & Middelberg targeted all kinds of FIRST sources with¬ 


out any preselection at a varying detection sensitivity. In order 


to compare their detection fraction with that of Garrtngton et al. 


Deller & Middelberg cut their catalogue to FIRST sources with 
> lOmJybeam * and to the VLBI detection sensitivity 


' 1.4 GHz 


from Garrington et al. Deller & Middelberg found a detection 
fraction of 36% in that subsample, in agreement with the num¬ 
ber from |Garrington et al.| 

We followed the approach from |Deller & Middelberg|p014) l 
and compiled a subsample of our IFRS sample by including only 
those sources with an arcsec-scale 1.4 GHz flux density above 
lOmJybeam ', ending up with 56 IFRS. Setting our detection 
sensitivity to 1.5 mJy, we would have detected 25 (45^g%) out of 
these 56 IFRS, i.e. a slightly higher fraction than those from Gar-| 


rington et al. (|1999 

1 andjDeller & Middelberg Using a Fisher's 

exact test (e.g. Wal 

& Jenkins 

2012 |l, we found a probability of 

0.88 that our samp 

e has a hig 

ler VLBI detection fraction than 


the sample from|Deller & Middelberg 


It is known that the radio source population at 1.4 GHz 
with flux densities above 1 mJy consists almost exclusively of 
AGNs (e.g. Condon et al.||2012| Fig. 11). This implies that the 
VLBI-observed subsamples from Garrington et al. ( 1999| l and 
[Deller & Middelberg] ( |2014| l and the sample presented in this 
work—all cut to lOmJy and matched to the same sensitivity as 
discussed above—contain virtually only AGNs. Thus, based on 
the numbers given above, we find a tendency of a higher VLBI 
detection fraction for IFRS compared to the general AGN popu¬ 
lation. 

[Collier et al.j ( |2014| l discarded all objects from their IFRS 
catalogue which provided a spurious WISE counterpart to one of 
the radio lobes. It is unclear whether this selection criterion can 
explain the higher VLBI detection fraction of IFRS compared to 
the general AGN population. Apart from that, the only difference 
in selecting the objects of the general AGN sample and the IFRS 
sample is the application of the IFRS selection criteria from [Zinn| 
[et al.| ( |20TT|l mentioned in Sect. |Ll| 

Figure [^ shows the VLBA detection fraction binned in the 
3 A pm flux density for our sample of IFRS. We do not find any 
evidence for a dependence of the VLBI detection fraction on the 
3 A pm flux density. We also tested the detection fraction against 
the arcsec-scale 1.4 GHz flux density, the radio-to-IR flux den¬ 
sity ratio, and the WISE colour [W1-W2] and found no signifi¬ 
cant correlation. 

Compact radio cores are detected in eight (73j;jg%) of the 


eleven lERS with SDSS counterpart and in 24 (57j^g%) of the 
42 lERS covered by SDSS without SDSS counterpart. However, 
we do not consider this a significant difference (a) because these 
two subsamples are not flux-complete, (b) because of the non- 
uniform sensitivity of SDSS, and (c) because SDSS is biased 
towards detecting blue objects as discussed above. 

Three out of the eleven lERS with SDSS counterpart are clas¬ 
sified as “galaxy” in SDSS. We detected two (67^2^%) of these 
three galaxy-type lERS in our VLBA observations. The other 
eight lERS with SDSS counterpart are classified as “star” in 
SDSS. We note that this photometry-based classification is based 
on the extension of the object, i.e. objects classified as star are 
point-like, whereas extended objects are classified as galaxy. Out 
of these eight lERS classified as star, we detected six (75j;jg%) 
in our VLBA observations. [Deller & Middelberg[(|20T4)l found a 
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Fig. 4. Histogram of the VLBA detection and non-detection fraction, 
binned by the flux density at 3.4pm. The horizontal lines show the over¬ 
all detection and non-detection fraction of the observations presented in 
this work. Icr confidence intervals are shown in grey. 


higher VLBI detection fraction for sources classified as star-like 
in SDSS. Our results are in agreement with this finding. 


4.1.2. Dependence of the VLBI detection fraction on radio 
properties 

Gigahertz peaked spectrum (GPS) and compact steep spec¬ 
trum (CSS) sources are very compact AGNs and expected to be 
the earliest phases in the evolution of AGNs. GPS sources have 
a turnover frequency of around 1 GHz and are usually less than 
1 kpc in size, whereas CSS sources are more extended with a 
size of a few kpc or a few tens of kpc. CSS sources are named 


for their steep radio spectra {a < -0.5; e.g. Randall et al. 


Based on data at 6 cm, 20 cm, and 92 cm. Collier et al. 


pOT^ 


classified 124 of their lERS as CSS sources and 32 as GPS 
sources. Out of the 57 lERS observed with the VLBA, five lERS 
were classified as CSS sources and two as GPS sources. In our 
VLBA observations, we detected four out of five lERS which 
were classified as CSS sources and both lERS which were clas¬ 
sified as GPS sources so that CSS/GPS sources have a higher 
detection rate (85^^^^ than the non-classified sources {5%^%). 

[Collier et al.|P()14] l used the lower-resolution data of NVSS 
for the flux densities in their lERS catalogue. They also listed 
the number of sources in the higher-resolution EIRST survey as¬ 
sociated with the NVSS source. lERS detected with more than 
one EIRST component are clearly extended radio galaxies and 
not GPS or CSS sources. Out of the 57 lERS observed with 
the VLBA, 47 are associated with exactly one EIRST source, 
while ten lERS are associated with two or three EIRST sources. 
Out of these ten lERS with two or three EIRST counterparts, 
we detected one in our VLBA observations, whereas 

34 (72+2*^) out of 47 lERS with exactly one EIRST countei-part 
were detected with the VLBA. We found a statistically signifi¬ 
cantly higher VLBA detection fraction for lERS with exactly one 
EIRST counterpart compared to the detection fraction of lERS 
with more than one EIRST counterpart. 

We also compared our VLBA detection fraction of lERS with 
exactly one EIRST counterpart to the detection fraction found by 
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Fig. 5. Histogram of the VLBA detections and non-detections, binned in 
the compactness at 1.4 GHz. Compactness is defined as the ratio of mas- 
scale flux density to arcsec-scale flux density. The vertical lines show 
the mean compactness—from left to right—of all VLBA-undetected 
IFRS, of all VLBA-observed IFRS, and of all VLBA-detected IFRS, 
respectively. The former two lines represent upper limits as indicated 
by the horizontal arrows. 


[Deller & Middelberg| ( |2014) for the general radio source popula¬ 
tion. As described above, we matched the arcsec-scale radio flux 
density and the VLBA detection sensitivity to lOmJybeam ' 
and 1.5 mJy beam \ respectively. We found a detection fraction 
of 54;^^% for those IFRS with exactly one FIRST counterpart, 
compared to a detection fraction of 36% for the gen eral radio 
source population above lOmJybeam ' measured by Deller & 


Middelberg Thus, our sample of IFRS with exactly one FIRST 


counterpart is statistically different to the general radio source 
population, based on a probability of 0.01 in a Fisher’s exact 
test (e.g. Wall & Jenkins|[T012| l that the two samples are taken 
from the same parent population. 


4.2. Compactness 


We determined ratios of the integrated mas-scale flux density in 
the VLBA observations at 1.4 GHz to the integrated arcsec-scale 
flux density in NVSS at 1.4 GHz to fall between 0.86 and 0.014 
for the IFRS detected in our VLBA observations. On average, 
this ratio, which we refer to as compactness, was 0.33 + 0.23. 
Figure shows the number of detections binned by the com¬ 
pactness. 


The mean compactness of our detected IFRS of 0.33 + 0.23 
is lower than that of the two former VLBI detections of IFRS 
where [Norris et al.| ( |2007| l and [Middelberg et al.| ( [2008b[ ) de¬ 
tected 88% and 58%, respectively, of the arcsec-scale flux den¬ 
sity. This discrepancy may be due to small-number statistics or 
because our fluxes are measured on smaller scales than those 
from [Norris et al.[ and [Middelberg et al.j The restoring beam 


of the LBA observations presented by Middelberg et al. 

51.7 mas x 23.6 mas, i.e. 17 times larger compared to the me¬ 
dian beam of 14.3 mas x 4.7 mas in our VLBA observations. 


Norris et al. did not image their uv data because of the poor 


uv coverage, but their angular resolution was similar to that of 



Fig. 6. Compactness as a function of redshift for those IFRS with red- 
shift information. VLBA-detected IFRS are shown by crosses and the 
6.15a upper limits on the compactness of VLBA-undetected IFRS are 
shown by arrows. Three spectroscopic (red markers) and seven photo¬ 
metric (black markers) redshifts were used in this analysis. Note that 
the errors on the compactness are around 20%. 


Middelberg et al.| Therefore, a lower fraction of detected flux in 
our VLBA observations could be expected. 


4.2.1. Dependence of the compactness on the redshift 

In the following, we tested our data against a potential corre¬ 
lation between redshift and compactness. SDSS DRIO provides 
spectroscopic redshifts for three out of our 57 IFRS. Two (z = 
2.11 and z = 2.62) of those are undetected and one (z = 2.55) 
is detected in the VLBA observations. However, IFRS F0912 at 
z = 2.62 was observed for only ~ 25 min with the VLBA, re¬ 
sulting in a sensitivity only half that of the other sources. Since 
the number of objects in this subsample is too low to test our 
data, we extended our subsample by including those IFRS with 
photometric redshifts presented in Sect. |3.2[ 

Figure|^shows the compactness as a function of the redshift 
for all ten VLBA-observed IFRS with redshift information. The 
data do not provide compact objects at low redshifts, whereas 
compact objects were found at higher redshifts. The data sug¬ 
gest a possible correlation between compactness and redshift 
for IFRS. We tested the data using a Spearman rank correla¬ 
tion test (e.g. [Wall & Jenkins[|2012| l and found a correlation 
coefficient between 0.66 and 0.52, indicating a positive corre¬ 
lation between redshift and compactness. A correlation coeffi¬ 
cient of - 1-1 and -1 represents an ideal correlation and anticor¬ 
relation, respectively, whereas an uncorrelated data set is repre¬ 
sented by a coefficient of 0. In our case, the probability that the 
two parameters are uncorrelated is between 0.019 and 0.063. The 
margin arises from the unknown compactnesses of the VLBA- 
undetected sources for which only upper limits are known. We 
determined this margin using a permutation test. Based on the 
strong positive correlation coefficients, we suggest a correlation 
between compactness and redshift for our sample of IFRS. When 
considering only the VLBA-detected sources, we found a corre¬ 
lation coefficient of 0.96 and a probability of 4.5 x 10 that the 
parameters are uncorrelated. We cautiously note that the putative 
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Fig. 7. Compactness as a function of the arcsec-scale 1.4 GHz flux den¬ 
sity. VLBA-detected IFRS are shown by crosses and the 6.75cr up¬ 
per limits on the compactness of VLBA-undetected IFRS are shown 
by black arrows. The grey horizontal lines represent—from top to 
bottom—the mean compactness of all VLBA-detected IFRS, of all 
VLBA-observed IFRS, and all VLBA-undetected IFRS, respectively. 
The latter ones represent upper limits as indicated by grey arrows. The 
black line shows the minimal detectable compactness depending on 
the arcsec-scale flux density, based on an mean detection sensitivity of 
450yuJy in our VLBA observations. 


positive correlation seems to be mainly based on the two highest- 
redshift IFRS in Fig. which might be outliers. Therefore, we 
emphasise that this suggested correlation needs further testing. 


4.2.2. Dependence of the compactness on the 1.4 GHz flux 
density 


Figurej^shows the compactness as a function of the arcsec-scale 
1.4 GHz flux density and includes detections and upper limits for 
the non-detections. We did not And compact radio-bright IFRS, 
whereas compact radio-faint IFRS are common in our sample. If 
we divide our sample at an arcsec-scale flux density of 60mJy, 
we And twelve (25^^^%) sources with compactnesses above 0.4 
and 36 sources with compactnesses below 0.4 in the 

fainter subsample. At arcsec-scale flux densities above 60mJy, 
we And nine (100_i7%) sources with compactnesses below 0.4 
and no (0^'^%) source with a compactness above 0.4. This is in 
agreement with results from Deller & Middelberg ( |2014 1 , who 
found a statistically significant anti-correlation between com¬ 
pactness and arcsec-scale 1.4 GHz flux density in their sample 
of randomly selected radio sources. We used a Spearman rank 
correlation test (e.g. [Wall & Jenkinsp012| l to test for putative 
correlations. However, because of the significant fraction of up¬ 
per limits in this plot, we can only narrow down the correlation 
coefficient to a rather broad range. We found that the correla¬ 
tion coefficient is between 0.39 and -0.15. Based on this test, 
we cannot exclude either a positive or a negative correlation or a 
decorrelation. 
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Fig. 8. VLBA map of IFRS F0030. The brighter component (upper left) 
has a complex morphology and a flux density of 21.8 mJy. The sec¬ 
ond component (lower right) has a flux density of 6.0 mJy and is sepa¬ 
rated by 442. mas, corresponding to a distance between 1.7 kpc and 
3.8 kpc at any reasonable redshift. 


4.3. Individual sources 

In the following, we discuss three individual sources which are 
of particular interest. 


4.3.1. F0398 


The only VLBA-detected IFRS with spectroscopic redshift is 
F0398 at z = 2.55, showing an arcsec-scale 1.4 GHz integrated 
flux density of 62.9 mJy. This corresponds to a /f-corrected 
1.4 GHz rest-frame luminosity of 2.3 x 10^^ WHz ', using the 
radio spectral index a - -0.72 between 20 cm and 92 cm from 
ICollier et al.| (|2014| l. In our VLBA observations, the source— 
shown in Fig. [IJ —is slightly resolved with a peak flux density of 
13.7 mJybeam^' and an integrated flux density of 17.2 mJy, cor¬ 
responding to a luminosity of 6.3x10^® W Hz ' on scales smaller 
than 146 pc X43 pc. Based on this luminosity, F0398 can be clas¬ 
sified as Fanaroff-Riley (FR; Fanaroff & Riley|1974| l type 11. The 
source has a compactness of 0.283. 


4.3.2. F0030 


A particularly interesting source is F0030 which has two spa¬ 
tially separated components in the VLBA map shown in Fig. 
The first, brighter component shows a mas-scale flux density of 
21.8 mJy (S/N = 55) and is spatially resolved with a complex 
morphology which is unique in our observations. The second 
component shows a flux density of 6.0 mJy (S/N = 46) and 
is separated by 442.1/^y3 mas. The linear distance between both 
components is between 1.7 kpc and 3.8 kpc at any redshift in the 
range 0.5 < z < 12. 


In order to obtain a rough spectral index of the components, 
we separately imaged the four lower-frequency basebands and 
the four higher-frequency basebands centred at 1.380 GHz and 
1.508 GHz, respectively. Fluxes were measured as described in 
Sect. 3.1 We obtained spectral indices of -1.2 + 1.2 and -1.2 + 
0.6 for component 1 and 2, respectively. 
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4.3.4. Could double component sources be a coincidence? 

Sources F0030 and F0257 have both two compact components 
separated by 445 mas and 155 mas, respectively. In the follow¬ 
ing, we estimate the probability that two close by components 
are unrelated background sources. NVSS found 1.8 million ra¬ 
dio sources in an area of 3 x 10"^ deg^, out of which less than 50% 


Fig. 9. VLBA map of IFRS F0257. The brighter component (left) has 
a flux density of 5.7 mjy. The second component (right) has a flux den¬ 
sity of 2. ImJy and is separated hy 155.51^g2 mas, corresponding to a 
distance between 0.6 kpc and 1.3 kpc at any reasonable redshift. 


4.3.3. F0257 

Source F0257—shown in Fig. |9] —consists of two individual 
components. While the brighter component shows a flux density 
of 5.7mJy (S/N = 77), the weaker component is 155.5 );q 2 mas 
distant and shows less than half of the other component’s flux 
density (2.1 mJybeam^', S/N = 36). This angular distance cor¬ 
responds to a linear distance between 0.6 kpc and 1.3 kpc. We 
measured mas-scale radio spectral indices between 1.380 GHz 
and 1.508 GHz of -0.3 + 0.4 and 0.2 + 0.9, respectively, for com¬ 
ponent 1 and 2. [Collier et al. (20141 measured an arcsec-scale 
radio spectral index of -0.6 between 6 cm, 20 cm, and 92 cm. 


We also tested our data for different potential correlations 
in Sect. [4 Although not all of them are statistically significant, 
in the following, we explore the astrophysical consequences if 
these results are confirmed by subsequent observations and de¬ 
scribe how these hypotheses can be tested. 

In Sect. 4.2. 2[ we found a tendency that radio-brighter IFRS 
are less compact. Deller & Middelberg ( 2014| l found the same 
behaviour when testing the general AGN population and argued 
that this anti-correlation might be explained by Doppler boosting 
effects as presented by |Mullin et al.| ( 2008| l. Mullin et al. studied 
a complete sample of narrow-line and broad-line radio galaxies 
and found an anti-correlation between radio luminosity and core 
prominence. They argued that higher-luminosity sources have 
higher boosting factors, associated with narrower boosting solid 
angles and a higher fraction of sources for which a Doppler sup¬ 
pressed core is seen. In contrast, lower-luminosity sources have 
lower boosting factors and wider solid angles, corresponding to 
a lower core supression fraction and a higher compactness. 

We argue that another factor might contribute to the observed 
behaviour in our sample. It is known that AGNs evolve from 
very compact to extended objects by forming and expanding 
jets, associated with an increasing total luminosity. Namely, GPS 
sources are most compact and evolve into CSS sources and fi¬ 
nally into the largest radio galaxies, FR type I or II (e.g. |0’D^ 
1998[ Snellen et al.|1999[ Fanti|2009] ). Following this sequence. 


we would expect younger AGNs to be more compact than old 
AGNs. 

We note that both effects might overlap since they both pre¬ 
dict a lower compactness at higher luminosities. Based on the 
slightly higher VLBI detection fraction of IFRS—which is ex¬ 
pected to be a result of higher compactnesses—compared to the 
general AGN population as found in Sections 4.1.1 [ and 4.1.2| we 
suggest that IFRS are on average younger than the general AGN 
population. This would be in agreement with results by |Midd^ 
berg et al. (201 lb]l and Collier et al. (|2014|l who found some 


have compact cores detectable in VLBI observations (Deller 
|& Middelberg|[2()14| l. This corresponds to a sky density of 
1.8 X 10"° arcsec The probability of finding an additional un¬ 
related source at a given position in an area of 0.5 " x 0.5 " is 
therefore of the order of 5 x 10“^. Thus, we can effectively rule 
out any chance that the two components found both in F0030 
and F0257 are physically unrelated. 


5. Discussion 


IFRS to be GPS and CSS sources. The higher VLBI detection 
fraction for IFRS classified as GPS and CSS sources compared 
to non-classifled IFRS reported in Sect. [4.1.2 also agrees with 
this reasoning. 

We suggest that—at least some—IFRS are young AGNs. 
However, the sample presented by [Collier et al.|p014j l contains 
IFRS with different characteristics. IFRS, which are associated 
with more than one FIRST source, are clearly no GPS or CSS 
sources but extended radio galaxies. Those sources would be 
expected to be older and less compact than IFRS with exactly 
one FIRST counterpart. We found evidence for this expected be¬ 
haviour in Sect. 14.1.21 based on a lower detection fraction for 
IFRS that are extended in arcsec-resolution images. We there¬ 
fore suggest that these extended IFRS are on average older and 
more evolved than the VLBA-detected IFRS. 

We found lower mean compactnesses of our VLBA-detected 
IFRS compared to the other two VLBI-detections of IFRS by 
Norris et al. ( 2007[ ) and Middelberg et al. ([2008b[ l. In addition 


Our observations increase the number of VLBI-detected IFRS 
from 2 to 37. Based on our detection fraction of 61^j% and 
a reasonably large sample size, we And strong evidence that 
most—if not all—IFRS contain AGNs. This result confirms ear¬ 
lier studies b y e.g. [Garn & Alexander[ ( [2008[ l, [Middelberg et al.[ 
( |201 lb[ ), and Herzog et al. ( [2014[ l, suggested compact cores 
in IFRS based on SED modelling, radio-to-IR flux density ratios, 
and emission lines in optical spectra. With higher sensitivity, we 
would have most likely detected more sources in our VLBA ob¬ 
servations. 


to the technical explanations given in Sect. [4.2[ this discrepancy 


can also be explained by boosting effects. The IFRS from Norris 


et al. and Middelberg et al.[ are radio-fainter than the IFRS anal¬ 
ysed in our work. Following the reasoning by [Deller & Middel- 


berg p014|l that brighter objects are more likely to be Doppler 


suppressed, lower compactnesses for our IFRS compared to the 
fainter IFRS from [Norris et al.[ and [Middelberg et al.[ could be 
expected. 


In Sect. 4.2.1 we found a statistically significant correlation 


between redshift and compactness, with higher-redshift IFRS 
being more compact. Two arguments can explain this correla- 


Article number, page 9 of 13 














































































































A&A proofs: manuscript no. document 


tion. (a) Higher-redshift objects have a tendency to be younger 
than low-redshift versions of the same class of object. Com¬ 
bining the increasing luminosity of GPS and CSS sources with 
time and the boosting-related argument of decreasing compact¬ 
ness with luminosity, higher-redshift—and thus younger and 
fainter—IFRS would be expected to be more compact, (b) At 
higher redshifts, IFRS are more likely to be located in gas-rich 
environments as shown for high-redshift galaxies (Klamer et al. 
2006 |l. The higher gas density confines these objects and keeps 


them more compact. 

These results are all in agreement with the scenario that IFRS 
are younger and therefore less luminous compared to the general 
AGN population, resulting in higher compactnesses and higher 
detection fractions. However, we stress that this putative con¬ 
nection between the age of IFRS and their VLBI properties is 
not statistically significant and needs further testing. 

Alternatively, the slightly higher VLBI detection fraction of 
IFRS compared to the general AGN population could also be 
explained by a higher dust content of IFRS, making it harder 
for the jets to expand and resulting in a more compact object. 
However, no evidence has been found that IFRS are obscured 
by dust. On the contrary. Collier et al. P014| l and Herzog et al. 
( |2014| l argued that the IR faintness of IFRS is not caused by dust 
extinction. The SEDs resulting from our photometric redshift fit¬ 
ting presented in Sect. 3.2 also indicate that at least some IFRS 
are very blue and do not support the hypothesis that a significant 
fraction of IFRS is associated with dusty galaxies. 

Future observations will help to test the hypotheses made 
in this work. In particular, additional VLBA observations— 
similar to the observations presented in this work—scheduled for 
semester 15A will increase the sample size of VLBI-observed 
IFRS, providing the basis for more robust tests. 

We plan to match arcsec-resolution radio data at higher and 
lower frequencies, enabling the measurement of spectral indices 
and turnover frequencies of IFRS. This information will bring 
out a putative overlap between IFRS and GPS/CSS sources and 
provide further insights into the evolutionary status of IFRS 
(Herzog et al., in prep.). 

Radio observations—exposing the intermediate-resolution 
morphologies of IFRS—could discriminate between the two 
mechanisms which may be responsible for the lower compact¬ 
ness of brighter objects: beaming and age. If beaming is the 
predominant cause for this effect, the radio-fainter IFRS should 
mainly be one-sided objects (core-jet) since a beaming effect is 
more likely to be seen for fainter objects as discussed above. 
In contrast, radio-brighter IFRS would be expected to show a 
more symmetric structure. However, if the anti-correlation be¬ 
tween compactness and arcsec-scale radio flux density is mainly 
driven by the age of the objects, no difference in the morphol¬ 
ogy would be expected, although radio-brighter objects should 
be larger. 

Two IFRS stand out from our sample since they were found 
to be composed of two mas-scale components, separated by a 
few hundred mas. There are four different explanations for those 
sources which we now discuss. 

(i) GPS/CSS double lobe sources: It is known that GPS and 
CSS double lobe sources can appear as separate components in 
VLBI observations ( [Snellen et al.||20()3] l. In that case, the two 
components would be hot-spots in the two jets and steep spec¬ 
tral indices would be expected (e.g. Hovatta et al.|2014] l, while 
the VLBI-undetected core of the source would be between both 


components. GPS sources are usually smaller than 1 kpc (O’Dea 
[1998^ , whereas CSS sources show extensions of a few to a few 
tens of kpc ( |Fanti|2009l [Randall et al.|2011| l. 
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(ii) Compact core and jet of a GPS or CSS source: Related to 
the first scenario, the two components could be the compact core 
of the AGN and a hot-spot in one jet. In that case, the spectral 
indices of both components can be different. While the hot-spot 
should provide a steep spectrum as discussed above, the core 
component would most likely provide a flat spectrum. However, 
the core spectrum could also be steep ( [Hovatta et al.|2014) l. 

(iii) Gravitational lensing: The appearence of two compo¬ 


nents can also be explained by gravitational lensing (e.g. Porcas 
1998|l. In that case, the emission seen as two components would 


originate from one distant source whose emission is deflected by 
the gravitation of a nearby object. Therefore, similar spectral in¬ 
dices for the two components are expected. In F0030 and F0257, 
the two components are too close to find a potential gravitational 
lensing effect in optical images. 

(iv) Binary black hole: The two components could also be 
a binary black hole (e.g. [Burke-Spolaor|20TT] ). In that case, the 
spectral indices of the components could be flat or steep or mixed 
as discussed in (ii) for the compact core. 

Based on the available data, we are not able to exclude any of 
those four different explanations because of the large error bars 
on the spectral indices. F0030 is unlikely to be a GPS source 
because of the linear size of more than 1.7 kpc. Following the 
correlation between intrinsic peak frequency and linear size of 
compact AGNs presented by [O’Dea & Baum[ \\991) , F0030 
would be expected to show a rest-frame turn-over frequency of 
less than a few hundred MHz. Particularly, F0030 has different 
characteristics than the high-redshift (z = 5.774) steep spectrum 
source J0836-I-0054 ( [Petrie et al.[[2003] Frey et ar][2005| l. This 
RL quasar shows arcsec-scale properties (5 i.4ghz = 1.75 mJy, 
a - -0.8) similar to those of IFRS. However, in contrast to 
F0030, J0836-I-0054 has a second arcsec-scale radio component 
which was undetected in VLBI observations and is most likely 
associated with a lower-redshift galaxy. 


6. Conclusion 

We observed 57 IFRS with the VLBA and detected compact 
emission in 35 of them. Based on these observations, we draw 
the following conclusions. 

- We tested the hypothesis that IFRS contain AGNs. Our ob¬ 
servations finally confirm the suggested compact cores in 
the majority of—if not all—IFRS, establishing IFRS as a 
new class of AGN. Our data increase the number of VLBI- 
detected IFRS from 2 to 37. 

- Our data suggest that radio-brighter IFRS are on average less 
compact. This finding agrees with the evolutionary scenario 
that young AGNs evolve by expanding jets, becoming radio 
brighter and less compact with time. However, boosting ef¬ 
fects may play a role, too. 

- We found a marginal tendency for IFRS to show a higher 
VLBI detection fraction compared to randomly selected 
sources with mJy arcsec-scale flux densities, i.e. mainly 
AGNs. In our sample, the detection fraction is higher for 
IFRS with exactly one FIRST counterpart and for IFRS clas¬ 
sified as GPS and CSS sources. 

- A statistically significant correlation between redshift and 
compactness was found in our data for IFRS with higher- 
redshift sources to be more compact. This is in agreement 
with higher-redshift sources being located in denser environ¬ 
ments and having a tendency to be younger. 

- Two sources show two components each, separated by be¬ 
tween 0.4 kpc and 3 kpc at any reasonable redshift. These 
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components might be jet/jet or core/jet components of an 
AGN, a binary black hole, or arising from gravitational lens- 
ing. 

All our findings are in agreement with the scenario that IFRS 
contain young AGNs which are in an early stage of their evolu¬ 
tion. Their jets are not yet formed or expanded significantly, re¬ 
sulting in a very compact source. When evolving, the jets expand 
and the total radio fluxes of the sources increase, while the com¬ 
pactnesses decrease at the same time. We note that some IFRS 
already formed jets as known from arcsec-resolution maps. 

Our analyses in this work were limited because of the low 
number of objects in relevant subsamples. Based on new data 
from VLBA observations of IFRS in semester 15 A and planned 
lower-resolution observations, we are aiming at extending our 
study and further testing the hypotheses presented in this work. 
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Table 1. Component catalogue of 57 IFRS observed with the VLBA. Listed is the identifier, the NVSS ID, the FIRST position, the NVSS peak and integrated flux densities, the VLBA peak 
flux density and integrated flux density, the VLBA S/N, and the VLBA beam size. If a source is found to be composed of several components in the VLBA observations, we list the components 
individually. In case of a non-detection, we quote a 6.75a- peak flux density upper limit. The identifier follows the numbering by|Collier et al.|j2014J. NVSS and FIRST data are taken from|Collier| 
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